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0 .  Introduction 


Many  problems  in  continuum  mechanics  involve  an  incompressi¬ 
bility  condition,  usually  in  the  form  of  a  divergence  constraint. 
The  numerical  discretization  of  such  a  constraint  presents  some 
interesting  problems  with  regard  to  stability.  As  an  important 
example  we  consider  the  two-dimensional  Stokes  equations 

-AU  +  VP  =  F  in  0  £  (R 2 , 

(1) 

7-U  =0  in  Q, 


with  appropriate  boundary  conditions  on  dQ .  This  has  the 
standard  weak  formulation 


Find  U  e  V  s  [H1 (Q ) ] 2  and  P  e  V  c  L2(Q)  such  that 


(2) 


a(U ,v)  +  b(y ,P)  =  (F,v)  V  v  e  r 
b(U,q)  =  0  V  q  e  V. 


The  bilinear  forms  a  and  b  are  given  by 


a(U.v)  -  2 


Q 


y  e ij(y)e ij(H)  dx 

i.  j 


b(v , P )  = 


7 • v  P  dx , 


Q 


2  2 

and  (F,v)  denotes  the  usual  [L  (Q ) ]  inner  product.  The  tensor 


1  fd  d  1  For  J 

c  ,  ( y)  is  the  symmetric  derivative  Vj  +  ax- ViJ  '  The  8Paces  \ - 


ij 


Y  and  V  depend  on  the  boundary  conditions.  For  no-slip 

q  1 


boundary  conditions:  Y  =  (H1(0)]2  and  V  =  L2  (0 )  n  jj  q  =  0  j- , 

0  I  uy - 

for  stress-free  boundary  conditions:  Y  -  the  orthogonal 


•d 

1  C!l._ 


□ 

□ 


Distribution/ 
Availability  Codes 
^Avall  arxl/or 
Spec  led 


Dlst 


complement  of  {^.^(y)  *  0}  In  [H1(Q)]2  and  V  =  L2(Q).  A 
natural  discretization  of  (2)  consists  in  choosing  finite 
dimensional  spaces  Y^  e  Y ,  c  y  and  determining 

Ujj  e  Y ^  and  e  such  that 

(3)  a<UN,v)  +  b(v,PN)  =  (F,v)  V  v  e  Y^ 

b(UN,q)  =0  V  q  e  Vh* 


The  main  obstacle  in  connection  with  (3)  is  to  find  spaces  Y 


N 


and  VN  so  that  the  discretization  is  stable  and  at  the  same  time 
has  good  approximation  properties.  A  reasonable  requirement 
concerning  stability  seems  to  be 

(4)  HU  -  UNII  1  +  IIP  -  PNI|  2  *  cfmin  ||U  -  v||  j  +  min  IIP  -  q||  J  , 

H  L  ('yerN  H  qeV^  L  •' 

with  C  independent  of  the  dimension  variable  N.  It  is  well 
known  that  the  Babuska-Brezzi  condition 


(5) 


[  V • v  q  dx 

min  max  -  *  c  >  0, 

qevN\{o>  yerN  llyil  jllqll  2 

H  L 


with  c  independent  of  N,  is  sufficient  to  quarantee  (4)  (cf. 
[2],  [7]).  If  the  pressure  spaces  are  chosen  equal  to  V*rN> 

then  (5)  is  equivalent  to  the  requirement  that  the  divergence 
operator 

7  •  :  r  — ►  V 

N  N 

has  corresponding  right  inverses 


2  1 

that  are  uniformly  bounded  in  #(L  ;H  ).  In  this  case  (5)  is  both 


a  necessary  and  sufficient  condition  that  the  quasi-optimality 
estimate 

||U  -  UMii  ,  *  C  min  HU  -  v)|  . 

'  H1  v *rN  ~  ~  H1 

holds  for  arbitrary  admissible  force  P  (cf.  [13]). 

The  most  natural  low  degree  finite  element  spaces  often  fail 
to  satisfy  the  Babuska-Brezzi  condition,  except  on  very  special 
triangulations .  A  remedy  is  to  appropriately  enlarge  the  velocity 
space  or  to  deplete  the  pressure  space;  such  approaches  are 
analyzed  for  continuous  piecewise  linear  (bilinear)  velocities 
with  piecewise  constant  pressures  in  [6]  and  [11]  respectively. 

For  continuous  piecewise  quadratic  velocities  one  has  the  well 
known  Taylor-Hood  element,  with  continuous  piecewise  linear 
pressures  (for  the  analysis  leading  to  (5),  see  [4]  and  [14]). 
Enlarging  the  velocity  space  or  depleting  the  pressure  space  is 
also  in  general  necessary  for  cubic  velocities  and  quadratic 
pressures  (cf.  [13]). 

For  continuous  piecewise  polynomial  velocities  of  total 
degree  <  p,  p  >  4,  the  situation  is  quite  different.  For  an 
arbitrary  triangulation  the  range  of  the  divergence  operator 
acting  on  the  velocity  space  has  a  very  simple  characterization  — 
it  consists  of  all  piecewise  polynomials  of  total  degree  s  p  -  1, 
except  for  a  certain  constraint  at  so-called  singular  vertices 
(cf.  [12]).  Furthermore,  for  fixed  p  2  4,  the  divergence  operator 
possesses  maximal  right  inverses,  the  norms  of  which  are  bounded 
independently  of  the  mesh  size  h.  To  paraphrase:  the  condition 
(5)  is  satisfied  for  such  velocities  if  the  pressure  space,  V  , 


is  chosen  to  be  7 • f 


Using  the  analysis  in  [15]  we  are  able  to 


2  „1 


prove  that  the  same  right  inverses  have  S(L  ;H  )  norms,  which 


are  bounded  by  some  power  of  p,  for  fixed  h.  These  results 


(p  >  4)  are  valid  also  for  piecewise  polynomial  spaces  with  a 


homogeneous  Dirichlet  condition  on  dQ  (assuming,  of  course,  Q 


is  polygonal). 


In  this  note  we  show  with  a  few  examples,  theoretical  as  well 


as  computational,  that  it  is  not  in  general  possible  to  find 


maximal  right  inverses  for  the  divergence  operator,  acting  on 


entire  polynomials  of  degree  <  p,  the  norms  of  which  are  bounded 


2  1 

in  8(L  ; H  ),  uniformly  in  p.  We  discuss  both  spaces  of  total 


and  separate  degree  <  p,  as  well  as  spaces  with  and  without 


boundary  condi t ions . 


The  lack  of  uniformly  bounded  right  inverses  for  the  discrete 


case  Is  somewhat  surprising  when  compared  to  the  continuous  case: 


it  is  easy  to  see  that  there  exists  a  right  inverse  (7*)  which 


maps  Hs  =  7«(HS+1)2  boundedly  into  (Hs+1)2,  V  s  >  0.  A 


similar  result  holds  with  homogeneous  Dirichlet  boundary 


conditions,  even  for  non-smooth  (polygonal)  domains  Q  (cf.  [1]). 


Methods  that  use  high  degree  polynomials  to  approximate  the 


solution  to  the  Stokes  equations  are  quite  common,  whether  they  be 


variationally  based  spectral  methods,  or  collocation  based  pseudo- 


spectral  methods  (cf.  [10]).  Another  possibility  is  the  so-called 


p-version  of  the  finite  element  method  (cf.  [3]):  it  uses  a  rather 


coarse  mesh  (triangulation  or  lattice)  and  achieves  convergence  by 


including,  in  a  variational  formulation,  piecewise  polynomials  of 


high  degree  relative  to  this  mesh.  Even  though  the  Babuska-Brezzi 


t 


condition  may  only  be  satisfied  with  a  constant  approaching  zero  as 


some  negative  power  of  p,  these  methods  seem  to  have  (nearly) 


optimal  convergence  rates  as  far  as  the  velocity  is  concerned.  We 


briefly  return  to  an  explanation  of  this  (at  least  for  variational 


methods)  towards  the  very  end  of  this  paper.  In  the  special  case 


of  periodic  boundary  conditions  it  is  normal  to  consider  spectral 


or  pseudospectral  methods  based  on  trigonometric  polynomials  instead 


of  polynomials.  The  resulting  methods  are  much  more  likely  to  be 


uniformly  divergence  stable  (see,  e . g . ,  [8]),  however,  they  are 


restricted  in  their  applicability  due  to  the  boundary  conditions 


To  complete  the  introduction,  let  us  give  an  interpretation 


of  the  constant 


j  7  *  v  q  dx 

J  A 


min  max  - 

qe*N\{0)  ||  yH  jllqll  2 

H  L 


in  terms  of  the  associated  matrices,  when  V„  =  V-r„.  in  order  to 


do  so,  we  first  need  to  specify  our  choice  of  norm  on  [H  ]  :  of 


the  many  equivalent  norms  we  take 


l  *  Z  I  | 5x7-1  da  +  |J  2  da 


1  =  1  0 


Let  be  a  basis  for  and  let  ^^^=1  be  a  basls  for 


%  -  "V  The  matrices  B  =  .<=1  and 


C  =  ^Ckf^k=i  ,f=l  are  deflned  bV 


(6. a) 


ak€  vk'  1  <  f  -  N' 


(6.b) 


bk  1=  J  357  *k  *  357  U  d*  +  j  *k  ^  *  J  u  d*‘ 

i=»l  Q  GO 


1  s  k ,1  s  N, 


( 6 .  c ) 


ck€  =  [  d-'  1  s  4  M. 


A  is  the  discrete  representation  of  the  divergence  operator  and 


( B-  ,  -  )  and  (O,*)  represent  the  quadratic  forms  llvll  and 

H1 


|q||  _  respectively. 

Ei 


With  these  definitions  it  is  easy  to  see  that  is  the 


-1/2  t  1/2 

smallest  singular  value  of  the  N  *  M  matrix  B  AC  ,  and 


this  in  turn  is  the  square  root  of  the  smallest  eigenvalue  of  the 


positive  definite  symmetric  M  *  M  matrix 


D  =  C1/2AB_1ATC1/2, 


For  any  q  e  VN  let  (V*)N  q  *=•  denote  the  element  of  minimal 


H  -norm  that  has  q  for  its  divergence.  By  a  "worst  possible 


pressure"  (as  far  as  divergence  stability  is  concerned)  we  mean  a 


q_  e  V„,  ||q_j|  _  =  1 ,  on  which  the  "minimal  norm"  right  inverse 

U  JN  U  ^  £ 


M 

(V*)N  attains  its  operator  norm.  If  x  e  R  denotes  a  unit 


eigenvector  for  the  matrix  D,  corresponding  to  the  smallest 


eigenvalue,  \x„,  then 

N 


So  -  Z|c_1/2s)3  *)  e  *N 

j=l 


is  a  worst  possible  pressure. 


L 


1.  Theoretical  and  Computational  results. 


Let  R  denote  the  square  ( -1 , 1 )  *  (-1 , 1 ) .  We  first 

consider  polynomials  of  separate  degree  <  p,  1 . e . ,  the  velocity 
2 

space  is  (Qp)  where 

(9)  <2p  =  span{x“x”  :  0  s  m,n  <  p), 
and  the  corresponding  pressure  space  is 

(10)  7*(Cp)2  =  spanfx^x"  :  0  <  m,n  <  p,  m  +  n<2p). 

Note  that  we  use  p  as  a  subscript  instead  of  the  dimension 

2  -l 

variable  N  =  (p  +  1)  .  As  before  (7*  )  denotes  the  right 

P 

inverse  with  minimal  H1(R)  norm. 


Proposition  1 . 

-1  2  2 
The  operator  (7*  )  :  7-  (<2  )  — ►  (Q  )  ,  p  t  1,  considered 

P  P  P 

2 

as  an  operator  from  a  subspace  of  L  (R)  to  a  subspace  of 
(H1(R))2  satisfies 

cp  s  ll( 7  •  )~1H  s  Cp, 

p  8(L  ;IT) 

with  constants  0  <  c  and  C  independent  of  p. 


Before  giving  a  proof  of  Proposition  1  we  make  a  few 
observations  about  orthogonal  polynomials.  Let  denote  the 

Legendre  polynomial  of  degree  n,  with  the  standard  normalization 


I  i  <«  = 


It  is  not  difficult  to  see  that 


)  ‘n  *  2K7i'Wx>  -  Vl|x|»'  n  *  *• 


The  polynomial  £  may  be  written  as  a  telescoping  series 
[ n/2 ] -1 


X  (^n-2j  ^n-2(j  +  l))  +  {  /J n 

J  _ A  V 


even 


and  consequently 


[ n/2 ] -1 


<„<x,  =  J_  {  V23-1  +  {^:  n 


From  this  we  conclude  that 


1  [a^n] 


dx  is  of  the  order 


[n/2] -1 

5> 

J-o 


1  .e.  , 


18  of  the  order  n2 . 


Let  qn(x)  «  0  <  n  <  p  denote  the  polynomials 


qn(x)  =  €n(x) »  0  <  n  <  p 


rx 

.(xl  -  1  Vi 


(the  notation  should  properly  be  q^  ,  since  the  definition  of 
qn  depends  on  p,  but  we  drop  the  superscript  for  convenience) 
A  simple  computation  gives  that  the  polynomial 


r<x>  "  alv-2  +  3  f  Vl'  P  *  2' 


satisfies 


12  2  l 

2  _af_  +  , 

3|2P+1  ( 2p+l ) ( 2p-l ) 2 


|  ( r(x)  )  2  dx  <  6  +  - - - g 

L  P  ( 2p-3 ) ( 2p- 1 ) 2 


p—  3  2 

Since  Kn)£=1  u  {/  j)  are  mutually  orthogonal  in  L  ,  and 

rX 

since  they  are  also  orthogonal  to  ^p_2  and  J  ^p-i '  lt;  foli( 

-1 

that 


1  r  P 


^^ct^q^x}  dx  is  equivalent  to 


Z2  -12  -3 

a^Jn+l)  +  Op(p+l)  ,  with  constants 


that  are  independent  of  p. 


It  is  convenient  to  work  with 


{qm(xl)qn(x2)  :  0  '  ra,n  -  p'  m  +  n  <  2P> 


as  a  basis  for  V«(Q  )  •  Based  on  (15)  we  get  that  if 


*  Z  °»nWW’ 


m ,  n=0 
m+n<  2p 


then 


r  2 

(q)  dx  is  equivalent  to 


Z  +  Z<(an.p)2+(V)2,,"+irlp'3' 


with  constants  that  are  independent  of  p. 


we  have 


Given  q  e  7  •  (<2p) 


q  *  Z  “rniWW 


m,  n=0 
ra+n<2p 

for  some  set  of  coefficients  {a  } .  Define 

nn 


u  = 


Z  amn (J  ‘qm|qn(x2>  +  Z%>ljXlqm]  V*! 


rXl 

n<m<p  -1  m<p  -1 

rX2  v*  jX2 

-1  n<p  -1 


y  a  q  (x  )  [  q  +  q  (x  )  [  q 

Z—  mn  ^m  1  J  n  Z—  pn*p  1  J  n 


m<  n<p 


It  is  clear  that  u  e  (Qp)  *  with  7*u  =»  q.  It  remains  to 


estimate 

the  H-norm  of 

U. 

Using 

(16) 

we  immediately  get  that 

feUlll  2  + 

°Xl  1  L  (R) 

N§57u2»] 

2  5 
-  (R) 

CM  2  ; 

l/(R) 

also 

R 

«-« 

o 

aio 
-  3a 

00 

,  for 

P  >  2, 

and  therefore 


Concerning  5^~u1 : 


I  f  S  I  «  C||q|| 

R  L  <R> 


„a  „  2 

1(3  x-j  U1 II  2 

X2  1  L  i  R) 


*1 

n<m<p  -1 


,17)  '  |  Z°-n(f +Z°n.p[j\]Vl<X2’ 


*1 

4 

m<p  -1 


L2(R) 


m-l 


-3 


<  C  2 *  ( m+l ) 

m<p  n=0 


y  a  \U  }  +  a  «  , 

Z-.  mn  ^dx  n J  mp  p- 1 


L  (-1,1) . 


10 


because  of  the  identities  (11)  and  (12)  and  the 
2 

L  -orthogonality  of  the  Legendre  polynomials.  From  (17)  we  get 
by  means  of  the  triangle  inequality,  Schwarz  inequality,  and  the 
estimate  (13) 


Il|x7ull|22 

0X2  1  L  ( R ) 


<  C>  (m+1 )  ~3  >  (a  ) 2 


Z(amn)2*Zn2  +  (amp,2p_1  ' 


'ran'  Z— 
n=0  n=0 


the  right  hand  side  of  which  is  bounded  by 


Z(amn»2  * 


/  mn  /  . 

n<m<p  m<p 


Using  the  above  estimate  in  combination  with  (16)  we  get 


ll§x7ulM22  "  Cp2  Nq^2 

0X2  1  LZ(R)  L 


L2(R) 


The  same  estimate  holds  for  — uo II  o  •  In  summary  we  have 

0Xl  *  LZ(R) 

thus  established 

N!  j  s  Cp  ||  q  ||  2 

H1 ( R)  L  i R) 

and  since  (V*)  -1q  is  the  field  of  minimal  H1-norm  it  follows 
P 


that 


<**> II  2  1S  Cp ' 


To  verify  the  second  inequality  of  this  proposition,  consider 

q*  = 

for  some  fixed  polynomial  r  (of  degree  <  p  -  1).  As  a  basis 

for  Q.  we  choose 
P 

(qm(x  )q  (x  )  :  0  <  m,n  <  p) . 
m  1  n  £ 


ssa 


H 


For  an  arbitrary  u  e  (Q  )  there  exist  coefficients  {a  }  and 

—  p  nn 


{13 _ }  such  that 

inn 


U1  =  Z  “mnWW 


m,  n=0 


U2  = 


tr 

Z  "mnWW- 


m,  n=0 


d  d  * 

If  V*u  =  ^ — u  +  ^ — u  =  q  ,  then  we  must  necessarily  have 
(j  X  j  X  u  X2  2 


P  x 

Z  I  Vp  dX 


r1  r  2 

=  j  7*u(x2  ,x2)«p(x2)dx2  =  r(Xj)  j  -Cpdx. 


Due  to  (12)  and  (14)  it  follows  that 


r  P  -j 

Xampqm  =  (2p  "  1,r‘ 

'-m=0 


and  therefore 


U1  *  Zqn(x2>  2^Lan»n^m^xl  ^  +  9p(x2)  (2p-l)  (J  'r  +  c] 


After  differentiation  with  respect  to  x2 


- — u  =  Y  f~q  ]  (x  )  Y"' o  q  (x  )  +  i  , (x„ ) ( 2p-l 
axj  1  Z-*  [dx^nj  2  ran^m'  l'  p-1  2 


from  which  it  now  follows  (by  orthogonality)  that 


rXl 

)  (  j  r  +  c)  , 
-1 


(19) 


J  [dx7ulj  d*  *  2  ( 2p-l )  j  (J  r  +  c) 

R  -1  -1 


dx 


At  the  same  time 


!  2(2P-U  |  j  r  -  H  1  r  HI 

-1  -1-1  L  ( -1 , 1 ) 


.V2  =  2p7T  llrH22 

L^(R)  L  (-1,1) 


and  so  we  have  proven  that  for  any  u  e  with  V*u  =  q  one 

has 

Hall  i  *  cp  11**11  2 

H  (R)  L  ( R ) 

(for  fixed  r) .  This  verifies  the  lower  bound  on  the  norm  of 

(*• Jl1. 


We  now  consider  polynomials  of  total  degree  <  p.  Without 

2 

boundary  conditions  the  velocity  space  is  (P  )  ,  where 

Pp  =  span(x^X2  :  m  +  n  <  p), 

2 

and  the  corresponding  pressure  space  is  V*(f  )  *  ^p-l'  Firs* 

some  computational  results  for  the  domain  R  =  (—1,1)  x  (—1,1). 

As  a  basis  for  Pp  in  these  computations  we  pick  products  of 
integrals  of  Legendre  polynomials  (supplemented  by  the  constant 
function) : 


(20) 


1, 


-1 


y-2 

n 

-1 


0  <  m,n  s  p-1, 


Os  m,n  and  m+n  <  p-2. 


we  pick 


As  a  basis  for  V*(P  )  =  9 

P  p-1 

(21)  t  (x,K  (x„)  0  <  m,n  and  m  +  n  s  p-1. 

in  l  n  i 

The  top  plot  in  figure  1  shows  the  smallest  eigenvalue  of  the 

matrix  D  (as  defined  in  (7))  for  p  varying  between  1  and 

18.  The  eigenvalues  were  computed  using  two  EISPACK  subroutines: 

first  the  matrix  was  transformed  (by  orthogonal  similarity 

transformations)  into  a  tridiagonal  matrix  using  subroutine  TRED2 , 

then  the  eigenvalues  were  computed  by  the  QL  method  (an  obvious 

2  1 

variant  of  the  QR  method)  using  subroutine  TQL2.  The  8(L  ;H  ) 

norm  of  the  "minimal  norm"  right  inverse  is  the  reciprocal  square 

root  of  the  smallest  eigenvalue.  The  numbers  do  not  clearly 

indicate  whether  these  right  inverses  are  bounded  independently  of 

p  —  if  anything  they  seem  to  indicate  that  the  norms  grow  as 

p  — i  ®,  but  only  as  a  very  small  power  or  possibly  a  logarithm  of 

p  (the  corresponding  solid  line  was  computed  by  linear  regression 

—0  3975 

on  the  last  four  points,  it  is  proportional  to  p  ) .  We 

note  that  by  a  slight  change  in  the  proof  of  Proposition  1  one  can 

2  1 

show  that  the  present  right  inverses  are  bounded  in  8(L  ;H  )  by 
Cp,  however,  that  is  clearly  too  conservative.  The  situation  is 
more  clear  cut  in  the  case  when  homogeneous  Dirichlet  boundary 
conditions  are  prescribed.  The  velocity  space  is  then  (£  )  , 
where 


9  =  9  n  (vLn 

p  p  '  ‘dft 


0) 


and  the  pressure  space  is  V*(^p)  £  P  j .  A  simple  count  of 


dimensions  gives  that  7*(Pp)^  has  co-dimension  9  in  Pp_j, 


p  t  5.  V*(P  )  is  the  common  nullspace  for  the  following  nine 


linearly  Independent  functionals:  the  Integral  over  R,  point 

P-2 


fa  v 


evaluation  at  each  of  the  corners  of  dR  and 
at  the  center  of  each  side  of  dR.  As  a  basis  for  9 


evaluated 
we  take 


the  same  elements  as  in  (20)  except  for  those  in  the  first  line. 


and  those  corresponding  to  m  =  0  or  n  =  0  in  the  second  line. 
Instead  of  computing  the  matrices  A  and  C  by  using  a  basis  for 

7*  (^p)2  we  use  a  basis  for  n  J  q  =  o|-,  the  same  as  in 

(21)  except  for  the  constant  function.  The  only  effect  of  this 

in  terms  of  eigenvalues  and  eigenvectors  is  to  add  0  as  an 

eigenvalue  of  D  with  multiplicity  8,  p  £  5  (in  the  case 
o  2 

p  =»  4 ,  V»(j>p)  only  has  codimension  8  in  P  and  the 

corresponding  multiplicity  of  0  becomes  7).  The  lower  plot  in 

figure  1  shows  the  smallest  positive  eigenvalue  of  D  for  p 

varying  between  4  and  18;  the  eigenvalues  were  computed  as 

before,  using  EISPACK.  The  8 ( ; H* )  norm  of  the  "minimal  norm" 

—  1  °  2  °  2 

right  inverse  (7*)  :  V • (9  )  — >  (9  )  is  the  reciprocal  square 

P  P  P 

root  of  the  smallest  positive  eigenvalue.  The  numbers  clearly 

indicate  that  the  norms  are  not  bounded  uniformly  in  p.  Based  on 

these  numbers  it  is  reasonable  to  conjecture  that  the  norms  grow 

3/2 

at  least  linearly  in  p  and  at  most  as  Cp  (the  solid  line  on 

_2 

the  graph  of  eigenvalues  is  proportional  to  p  ,  corresponding 
to  linear  growth  of  the  norms).  What  is  more,  we  can  actually 
prove  the  lower  bound: 


I  r L4»  rv*  irwiTHTI 


Proposition  2 . 

Let  (V*)  *:  7*(^  )2  — ►  (^  ) 2 1  p  s  4,  denote  the  right 

^  r  r 

inverse  with  minimal  H*  norm.  Considered  as  an  operator  from  a 

2  12 
subspace  of  L  (R)  to  a  subspace  of  (H  (R) )  this  satisfies 


cp  s  II  (V-  )  "1|| 


B(L4 


H1) 


with  c  >  0;  independent  of  p. 


Proof . 


* 

Let  q  (x1(x2)  *  ^1(x1)  j  'ep_3-  I1:  is  clear 


that 


-1 
* 


fxi  rx2 

q  e  7  •  (^p)fc  since  ?  •  v  =  q  for  v  «=  ( j  ^  I 

-1  -1 


Vs-  0)-  0n 


the  other  hand,  if  u  e  (Pp)  is  an  arbitrary  velocity  field  with 
* 

V  •  u  =  q  ,  then 


u  = 


m+nsp-2 

y  a  \  2< 

Z-  an  J  m  J  n 
m , n=l  -1  -1 


m+nsp-2 

Zr  X1  fx2 

(i  *  \  t 

mn  J  m  J  n 

m , n=l  -1  -l 


and 


m+n<p-2 


I  (■ 


a  <  (x  ) 
mn  m  l 


.x2 


n 


+  fi 


mn 


Xi 


VW 


J" 


m , n=l  -1  -1 

Using  (12)  and  the  linear  independence  of  the 
mials  we  get 


-1 


Legendre  polyno- 


17 


vinu  wirwm  u 

Bf 


I 


& 


o 


g 


i 


P-3  x 

I«.n  J  «a 

n=l  -1 


Z  5  *2nW 


P-3' 


since  these  are  the  respective  coefficients  of  ^(Xj). 


identity  (22)  implies 


alp-3  =  lf 


and  since 


m+n<  p-2 


ilE7ul  =  X  amn(I  VW 


m ,  n=l 


m+n<  p-2 

^  2m+Iamn[*m+l  (xl  >  *m-l  (xl *]  €n(x2 }  ' 

m ,  n=l 


it  follows  that 


(23) 


»aV4,R)  *  lI(o  }% 

L  (R)  -1  -1 


9(2p-5)  * 


A  simple  computation  gives 


L  <R) 


8  1 
3  ( 2p-3 ) ( 2p-5 ) ( 2p-7 )  ' 


and  therefore,  in  light  of  (23), 


||u||2  *  |  (  2p-3 )  (  2p-7 )  || q4 


H  ( R) 


I  2 
L*(R) 


Since  u  is  an  arbitrary  field  in  (£p)2,  with  V  •  u 


this  gives  the  desired  lower  bound  on  the  operator  norm. 


Figures  2  and  3  show  contour  plots  of  worst  possible  pressures 

in  tf  ?  sense  defined  at  the  end  of  the  introduction,  for  p  =  7,  8, 

o  2 

14  and  15,  These  pressures  are  elements  of  V*(J»  )  ,  and  they 

have  the  property  that  the  right  inverses  (V*)  1:  7 •  (9  )2  — ►  (9  )2 

P  P  P 

2  1 

attain  their  ®(L  ;H  )  norms  there.  Solid  lines  in  the  plots 
correspond  to  positive  values,  dashed  lines  correspond  to  negative 
values.  On  each  plot  the  interval  between  contour  lines  is 
indicated  at  the  bottom  and  so  is  the  entire  value  range  of  the 
pressure . 

We  note  several  features: 

(1)  There  is  a  marked  difference  between  worst  pressures  for 
even  and  odd  p.  This  correlates  well  with  the  lowest  plot  in 
figure  1,  which  really  seems  to  consist  of  two  slightly  different 
curves,  one  for  even  p  and  one  for  odd  p. 

(2)  The  value  range  for  a  worst  pressure  grows  as  p  in¬ 
creases,  and  the  extreme  values  are  clearly  attained  on  dQ . 

(3)  For  p  odd  and  sufficiently  large  there  is  a  local 
checkerboard  pattern  developing,  similar  xo  that  found  in 
connection  with  some  unstable  low  order  elements  (see,  e . g . ,  [5] 
for  bilinear-constant  velocity-pressure  approximation) . 


2. 


Conclusions  and  final  remarks . 

We  have  shown  with  a  few  examples  that  one  cannot  In  general 


construct  maximal  right  inverses  for  the  divergence  operator, 

2  1 

whose  £(L  ; H  )  norms  are  uniformly  bounded  as  the  polynomial 
degree  increases.  We  do  not  think  that  the  square  domain  is 
extremely  special  in  this  regard,  and  we  believe  that  similar 
examples  can  be  found  on  many  other  polygonal  domains. 

We  do  not ,  however,  claim  that  it  is  impossible  to  find 
domains  for  which  uniformly  bounded  maximal  right  inverses  exist. 
Indeed  one  such  class  of  domains  { for  polynomials  of  total  degree 
<  p,  with  no  boundary  conditions)  are  the  ellipses: 

Example  2 . 1 . 

2  2 

Let  E  =  { ( x j , x^ ) :  aXj  +  bx2  <  1),  0  <  a,b.  The  Laplace 
operator  A  :  j^-j 2  +  j^~ -j 2  maps  fp+J  into  f  p-1 .  The  space 

o  2  2 

P  ^ .  ,  on  the  domain  E,  is  the  same  as  (ax.  +  bx_  -  1)9  ,, 

P+1  1  2  p-1 

and  since  A  has  no  nontrivial  null  vectors  with  homogeneous 
Dirichlet  boundary  conditions,  it  follows  that  A  is  an 
isomorphism  from  Pp+1  onto  Fp_1.  The  operator 
-l-l  2 

(7 • ) p  =  7 A  :  — ►  (Pp)  is  a  uniformly  bounded  maximal  right 

inverse  for  the  divergence  operator  on  the  domain  E. 

■ 

Where  no  boundary  conditions  are  involved,  unboundedness  of 
the  “minimal  norm"  right  inverse  in  spaces  of  entire  polynomials 
(on  a  square,  say)  immediately  leads  to  unboundedness  in  spaces  of 
piecewise  polynomials  relative  to  a  fixed  partition  (a  lattice). 

We  expect  that  the  "minimal  norm"  right  inverses  for  truly 


piecewise  polynomials  will  inherit  some  of  the  (likely)  extra 


unboundedness  associated  with  homogeneous  boundary  conditions 


For  example,  on  a  lattice  (with  more  than  one  rectangle)  it  would 


not  be  surprising,  if  the  right  inverses  corresponding  to  spaces 


2  „1 


of  piecewise  polynomials  of  separate  degree  <  p  have  8(L  ;H  ) 


norms,  that  grow  faster  than  p. 


Lack  of  divergence  stability  as  evidenced  by  the  fact  that 


the  best  lower  bound  in  (5)  might  behave  like  p  ,  a  >  0, 


surprisingly  does  not  lead  to  suboptimal  order  of  convergence  for 


the  velocities,  as  p  — >  ®  (provided  we  use  the  divergence  of  the 


velocity  space  as  the  pressure  space).  This  may  be  explained  by 


an  interpolation  argument,  the  idea  of  which  originates  in  [3). 


The  explanation  is  particularly  simple  with  stress-free  boundary 


conditions  on  a  simply  connected  domain:  since  Up  is  a  pro¬ 


jection  of  U  we  get 


1“  -  Spll  i 

P  H  (0) 


^  cllul 


H1(Q) 


and  since  the  lower  bound  in  (5)  is  -  p  we  can  also  prove 


(cf.  [13]) 


|U  -  U  ||  *  Cpa  min  ||U  -  y|| 

P  H  (Q)  v  H  (Q) 


where  the  minimum  is  taken  over  the  full  space  of  discrete 


2  2 

velocities  (  (<2  )  or  (?  )  )  .  Standard  approximation  results 
P  P 


yield 


IS  -  Spll  ,  «  *1 


,  Mi  0 , 


H  (Q  ) 


H  (Q  ) 


v'-,  .V  AUVV 


fMTt 


For  fixed  M  we  can  interpolate  a  fraction  0  =  n  between  (24) 
and  (25)  to  get 


la  -  a  ||  ,  c  p-k+ak/M  ||u|| 

P  H  (0)  H*  (0) 


,  0  <  k  <  M. 


By  choosing  M  sufficiently  large  (M  >  — ^)  we  conclude  that 


-k+£ 


=  -  VL 5  c*.£p 


H‘(Q) 


H  (Q  ) 


for  any  c  >0.  Modulo  £  the  estimate  (26)  represents  the 

lc+ 1 

optimal  order  of  convergence  for  general  U  in  H  (Q).  More 
details  are  found  in  [15]  and  [16],  where  the  same  argument  is 
used  in  the  context  of  the  equations  of  elasticity  (for  a  nearly 
incompressible  material). 

N 

We  do  believe  that  the  lack  of  divergence  stability  affects 
the  accuracy  of  the  pressure  approximation  much  more  drastically, 
and  we  expect  that  a  certain  postprocessing  (filtering)  of  the 
pressures  may  be  necessary  as  p  ® .  •  > 
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